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1.  Introduction 


The  U.S.  Army  Research  Laboratory  Gun  Tube  Erosion  Code  (ATEC)  is  the 
culmination  of  effort  by  the  U.S.  Army  Research  Laboratory  (ARL)  to  identify  the 
thermochemical  components  to  gun  tube  erosion.  The  inclusion  of  chemistry 
was  originally  conceived  due  to  discrepancies  of  experimentally  fired  propellants 
to  the  previous  correlations.  The  historical  correlations  were  based  solely  on  the 
propellant's  adiabatic  flame  temperature  and  did  not  take  into  account  the 
possible  chemistry.  To  address  this  deficiency,  ARL  embarked  on  a  project  to 
include  chemistry  in  a  model  which  would  enable  the  prediction  of  erosion  from 
basic  principles  without  requiring  too  much  information  prior  to  the  actual 
system  investigation.  The  model  needed  to  be  system  independent  with  enough 
basic  physics  such  that  it  would  function  for  most  guns  or  gun-like  systems 
irrespective  of  caliber.  To  this  end,  a  modular  gun  tube  erosion  code  was 
developed.  Technical  details  of  the  code,  as  well  as  results,  have  been  previously 
reported  by  Conroy  et  al.  [1-6]  and  Weinacht  and  Conroy  [7]. 

The  following  physical  phenomena  are  included  in  the  model: 

•  Convective  heat  transfer  to  the  surface  through  convection  coefficients, 

•  Subsurface  heat  conduction  with  variable  physical  properties  for  gun  steel, 

•  Multicomponent  species  mass  transport  to  the  surface,  independent  of  heat 
transport, 

•  User-specified  finite  rate  surface  reactions, 

•  Subsurface,  diffusion  limited,  equilibrium  reactions  with  user-specified 
reactants  and  products, 

•  User-specified  coating  materials, 

•  Coating  stress/ crack  distribution  calculation, 

•  Erosion  at  the  base  of  cracks  in  surface  coatings,  and 

•  Uncoated  surface  erosion. 

This  report  documents  the  ATEC  user  input  as  well  as  output  from  the  code. 
Instructions  for  developing  an  input  deck  are  explicitly  provided  in  the  following 
sections.  Files  required  to  operate  the  code  include  a  user-generated  input  deck, 
"erod.in",  an  output  file  "htemp.out"  from  XKTC  [8],  the  NASA  Lewis  library 
"library .in"  [9],  and  a  species  temporal  file,  "spedes.out",  from  a  coupled  interior 
ballistic-equilibrium  thermodynamic  code. 

Example  input  decks  are  provided  in  Appendices  A  and  B.  An  example  of  the 
"spedes.out"  file  is  given  in  Appendix  C.  The  ATEC  code  has  been  applied  to 
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the  25-mm  Bushmaster  cannon.  Navy  MK45  5-in  cannons,  the  120-mm  M256 
tank  cannon,  as  well  as  various  vented  fixtures,  and  has  performed  well  in 
representing  the  physics  occurring  at  the  surface. 


2.  General  Computer  Model  Description 


Figure  1  represents  what  occurs  beneath  the  coating  near  a  crack.  What  is 
evidenced  is  that  a  jet  of  combustion  gases  washes  through  the  crack 
subsequently  reacting  and  eroding  the  substrate  material  leaving  a  pit  with 
horizontal  extensions  near  the  interface.  Many  species  of  materials  in  the  eroded 
region  have  been  identified  and  include  sulfur,  oxygen,  and  talc  from  the 
propellant  additives  [10, 11].  An  interesting  observation  is  that  sulfur  penetrates 
deep  into  the  eroded  pits.  Sulfur  comes  from  the  black  powder  ignitors  and 
potassium  sulfide  flash  suppressants  used  in  the  charges. 

Figure  2  is  a  cartoon  of  the  numerical  simulation  of  erosion  at  the  bottom  of  a 
radial  crack.  It  is  important  to  realize  that  crack  orientation  is  critical,  for  if  the 
crack  is  longitudinal,  it  would  enable  propellant  product  gases  to  wash  directly 
in  and  out  with  tremendous  amounts  of  heat  transfer  and  the  possibility  for 
surface  reaction,  in  comparison  to  a  radial  crack  which  would  rely  more  on 
recirculation  and  diffusion  of  gases  to  the  interface. 

Figure  3  shows  the  interconnectivity  of  the  codes  for  an  erosion  calculation.  The 
interior  ballistics  is  modeled  yielding  the  thermochemistry  of  the  propellant 
products  for  further  transport  to  the  surface,  as  well  as  supplying  gas  state  and 
velocity  variables  for  convection  heat  transfer.  The  erosion  code  itself  has  an 
extensive  input  deck,  the  description  of  which  is  the  purpose  of  this  report. 


3.  Required  U.S.  Army  Research  Laboratory  Gun  Tube 
Erosion  Code  (ATEC)  Input  Data  Files 


Table  1  provides  a  guide  to  all  input/ output  files  for  the  ATEC  code.  The 
following  are  required  input  files: 


(1)  A  mixed  free  format  and  formatted  input  data  deck  specific  to  ATEC  called 
"erod.in".  This  input  deck  is  described  in  detail  in  this  guide. 

(2)  A  data  file  from  the  XKTC 1-D  NOVA  code  entitled  "htemp.out". 
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(3)  A  data  file  entitled  "species.out".  An  example  of  this  file  is  supplied  in  this 
report.  It  can  be  generated  by  any  number  of  coupled  interior  ballistic  and 
thermodynamic  equilibrium  codes  such  as  XKTC  with  NASA  Lewis. 

(4)  The  NASA  Lewis  species  thermochemistry  library.  This  provides  the 
thermochemistry  for  both  the  equilibrium  and  finite  rate  kinetics 
calculations. 

Note:  Lower  case  should  be  used  for  input  file  names. 

The  following  is  the  format  required  of  the  file  "species.out". 

The  first  line  contains  the  header  "-Time_(ms)"  followed  by  the  names  of  the  up 
to  top  25  resultant  species  from  the  equilibrium  calculation  entered  in  the 
following  format  where  the  names  are  entered  as  "6al"  in  standard  NASA  Lewis 
format. 

READ(30,850)  ((ZLABEL(II,JJ),JJ=1,6 ),  II=1,SPNDIM2) 

850  FORMAT ( '  -Time_(ms)',  25(',  -'MV-') ) 

Note:  The  species  should  be  in  descending  order  from  highest  initial  molar 
concentration  to  the  lowest,  from  the  second  to  the  leftmost  column  to  the  last 
column  to  the  right. 

Inclusion  of  the  top  25  species  should  suffice  for  calculations  purposes.  The  lines 
following  the  first  line  of  the  input  file  "species.out"  contain  the  time  history  of 
each  species  concentration  with  the  first  column  tracking  time  in  milliseconds 
and  the  following  columns  containing  the  species  concentration  in  mole  fractions 
in  the  format  presented  in  the  following  read  statement: 


READ(30,860,END=2000)(VAR(KT1,KN1,NV),NV=1,SPNDIM2+1) 
860  FORMAT(  1X,F10.4, 25(7,1PE10.3) ) 

An  example  of  a  "species.out"  file  is  provided  in  Appendix  C. 


4.  ATEC  Input  Deck  Generation  Guide 


0.1  Examine  the  htemp.out  output  file  from  XKTC. 

Determine  initial  grid  spacing  (in). 

Determine  initial  temperature  (R). 

Determine  initial  chamber  length  (in). 

Determine  final  length  of  the  tube  from  the  last  time  step  (in). 
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1.1  Enter  the  initial  temperature  from  the  htemp.out  file  (R). 

1.2  Enter  the  initial  pressure  from  the  htemp.out  file  (R). 

1 .3  Enter  the  initial  velocity  from  the  htemp.out  file  (in/s). 

2.1  If  there  is  an  inner  surface  coating  layer  enter  thickness  (in). 

2.2  Enter  total  time  of  each  firing  (ms). 

2.3  Usually  the  chamber  augmentation  factor  is  0.1  (-). 

2.4  Enter  the  chamber  length  (in). 

2.5  Enter  the  tube  length  from  the  last  time  step  of  htemp.out  file  (in). 

2.6  Enter  the  number  of  axial  grid  points  (<100). 

This  is  not  a  trivial  task! 

The  INITIAL  (t=0)  grid  spacing  horn  the  htemp.out  file  MUST  BE  GREATER 
than  the  erosion  grid  spacing.  If  you  wish  to  examine  the  entire  barrel  length, 
you  may  have  to  accomplish  this  in  more  than  one  calculation  as  the  code 
currently  accepts  100  axial  points,  where  NOVA  usually  has  50  or  more 
distributed  over  about  a  40-in  chamber,  for  example. 

2.7  Determine  the  computational  region  desired. 

Make  sure  the  length  of  the  computational  space  doesn't  expand  the  grid  spacing 
beyond  the  initial  htemp.out  grid  spacing. 

Enter  the  calculation  starting  position  (in). 

This  must  be  within  the  specified  tube  geometry. 

2.8  Enter  the  calculation  ending  position  (in). 

This  must  be  within  die  specified  tube  geometry. 

2.9  If  2-D  LP  chamber  calculation,  enter  a  one  to  trip  flag  (LP  option). 

(Not  currently  implemented  in  erosion  code,  i.e.,  must  be  zero.) 

2.10  If  2-D  LP  chamber  calculation,  enter  the  initial  piston  offset  (in) 
(LP  option).  (Not  currently  implemented  in  erosion  code,  i.e.,  must  be 
zero.) 

2.11  If  IMC  cooling  is  desired,  enter  1  (this  is  a  flag  for  an  extra  card). 

(Not  currently  implemented  in  erosion  code,  i.e.,  must  be  zero.) 

2.12  If  the  outside  surface  is  insulated,  enter  1. 

(Not  currently  implemented  in  erosion  code,  i.e.,  must  be  zero.) 

2.13  If  blowby  is  desired,  enter  1  (this  is  a  flag  for  an  extra  card). 

(Not  currently  implemented  in  erosion  code,  i.e.,  must  be  zero.) 
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2.14  Enter  the  amount  by  which  you  want  to  multiply  the  internal  convective 
heat  transfer  coefficient.  If  no  multiplication  is  desired,  enter  1.D0  or  0.0 
(both  produce  the  nominal  result). 

Cl.l  This  is  conditional  on  the  heat  transfer  coefficient  augmentation  in  2.14. 
If  a  negative  value  is  entered,  then  a  roughness  calculation  is  performed 
based  on  defects  on  tube  surface  and  Moody  chart.  Enter  smooth 
deformation  limit  (m). 

Cl. 2  Enter  rough  deformation  limit  (m). 

3.1  Enter  the  number  of  axial  locations  to  be  specified  (-). 

3.2-9  Enter  up  to  eight  axial  locations  (in). 

4.1- 8  Enter  up  to  eight  corresponding  outer  diameters  (in). 

5.1- 8  Enter  up  to  eight  corresponding  inner  diameters  (in). 

Breech  is  zero.  Calculation  occurs  in  specified  computational  space. 

6.1  Enter  the  number  of  probe  locations,  up  to  eight  axial  locations  (-). 

6.2- 9  Enter  up  to  eight  axial  locations  (in). 

7.1-8  Enter  up  to  eight  corresponding  probe  locations  in  the  tube  (in)  . 

Note:  The  probe  depth  is  only  printed  to  two  decimal  places,  so  don't  be 
surprised  when  you  ask  for  a  probe  at  0.005  inches  and  it  comes  back  as 
0.01  inches.  It  is  just  an  artifact  of  the  print  format.  The  probe  depth  is 
from  the  bore  surface  into  the  tube! 

8.1  Enter  "1"  to  activate  only  the  equilibrium  surface  chemistry. 

Enter  "2"  to  activate  both  the  subsurface  equilibrium  chemistry  and  finite 
rate  surface  chemistry 

Enter  zero  if  you  wish  to  bypass  the  chemistry  (-). 

8.2  Enter  the  freeze-out  temperature  (K),  i.e.,  temperature  at  which  the 
transport/  surface  chemistry  routine  is  not  called. 

8.3  Enter  density  of  the  base  wall  material  (kg/m3). 

8.4  Enter  temperature  at  which  the  base  wall  material  begins  to  melt  (K). 

8.5  Enter  specified  outer  barrel  heat  transfer  coefficient  (W / m2/K). 

8.6  Enter  ambient  outer  air  temperature  (K). 

8.7  Enter  the  computational  coordinate  spacing  (-)  usually  "1". 

8.8  Enter  the  latent  heat  of  fusion  for  the  material  (J/kg). 

8.9  Enter  "0"  if  the  system  is  planer  or  "1"  if  the  system  is  axisymetric. 
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8.10  Wall  temperature  calculation  flag  (-). 

Molecular  weighting  function  between  gas  and  solid. 

Wall  temperature. 

Conditional  Card  (if  Chrome  >  0). 

C2.1  Number  of  radial  points  considered  in  coating  for  calculation  (-). 

C2.2  Number  of  radial  points  considered  in  base  material  for  calculation  (-). 

C2.3  Enter  the  freeze-out  temperature  of  surface  material  (K),  i.e.,  temperature 
at  which  the  chemistry  routine  is  not  called. 

C2.4  Enter  density  of  the  coating  material  (kg/m3). 

C2.5  Enter  temperature  at  which  the  coating  material  begins  to  melt  (K). 

C2.6  Enter  the  computational  coordinate  spacing  (-)  usually  "1". 

C2.7  Enter  the  latent  heat  of  fusion  for  the  coating  material  (J/kg). 

C2.8  Stress  flag  (0)-no  stress  calculation  (l)-stress  calculation  (-). 

Conditional  Card  (if  Chrome  >  0). 

C3.1  Enter  the  Young's  modulus  of  the  coating  material  (Pa). 

C3.2  Enter  the  Coefficient  of  Thermal  Expansion  of  the  Coating  (1/K). 

C3.3  Enter  the  Poisson's  Ratio  Coating  (-). 

C3.4  Enter  the  ultimate  Strength  of  Coating  (Pa). 

C3.5  Enter  the  Young's  Modulus  Base  (Pa). 

C3.6  Enter  the  coefficient  of  Thermal  Expansion  Base  (1/K). 

C3.7  Enter  the  Poisson's  ratio  of  the  base  material  (-). 

Conditional  Card  (if  Chrome  >  0). 

C4.1  Enter  crack  width  (m). 

C4.2  Enter  crack  distribution  limit  (m). 

10.1  Enter  the  number  of  shots  to  fire. 

(Not  implemented  in  this  version  of  Erosion  Code,  i.e.,  must  be  1). 

10.2  Enter  print  frequency  during  htemp.out  file  period  (-).  Choosing  the 
computational  time  interval  (implicit)  (usually  about  0.1  ms)  enables  the 
print  frequency  to  be  entered,  i.e.,  print  every  x  time  steps. 

10.3  Enter  print  frequency  after  htemp.out  file  period  (-). 
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10.4  If  restarting,  enter  the  number  of  the  last  round  fired  before  restart,  else 
enter  a  zero. 

(Not  implemented  in  this  version  of  Erosion  Code,  i.e.,  must  be  0;. 


11.1-11.17  Control  the  files  to  be  printed. 

(o  =  false),  (1  =  true),  (2  =>single  shot  =  burst) 

1.  0  =  K,  1  =  R  (units  Rankine  or  Kelvin) 

2.  0  =  as  computed,  1  =  temperature  rise 

3.  Screen  output  IOUT(l) 

4.  Single  shot  probe  temperature  IOUT(2) 

5.  Single  shot  inside  surface  temperature  IOUT(3) 

6.  Single  shot  outside  surface  temperature  IOUT(4) 

7.  Burst  fire  probe  temperature  IOUT(5) 

8.  Burst  fire  inside  surface  temperature  IOUT(6) 

9.  Burst  fire  outside  surface  temperature  IOUT(7) 

10.  Heat  transfer  coefficient  IOUT(8) 

11.  Effective  gas  temperature  history  IOUT(9) 

12.  Single  shot  graphic  output  file  (for  IBGRAF  movies)  IOUT(IO) 

13.  Burst  fire  graphic  output  file  (for  IBGRAF  movies)  IOUT(ll) 

14.  Boundary  condition  graphic  output  (for  IBGRAF  movies)  IOUT(12) 

15.  Bore  surface  heat  flux  output  (single  shot-1,  burst-2)  IOUT(13) 

16.  Bore  surface  regression  rate  (single  shot-1,  burst-2)  IOUT(14) 

17.  Bore  surface  regression  (single  shot-1,  burst-2)  IOUT(15) 

18.  Bore  surface  concentrations  IOUT(16) 

19.  Bore  core  flow  concentrations  IOUT(17) 

20.  Shear  stress  azmuthal  at  interface  IOUT(18) 

21.  Shear  stress  axial  at  interface  IOUT(19) 


(Note:  Screen  output  is  for  PC  applications  only  and  is  not  supported.) 

12.1  Variable  4340  conductivity  and  diffusivity  flag  1-use  internal  functions, 
0-use  input  values. 

12.2  Barrel  conductivity  (W /m/K)  (Necessary  even  if  12.1  is  1.) 

12.3  Barrel  diffusivity  (m2/s)  (Necessary  even  if  12.1  is  1.) 
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12.4  Liner  conductivity  (W/m/K)  Necessary  even  if  12.1  is  1.) 

12.5  Liner  diffusivity  (m2/ s)  (Necessary  even  if  12.1  is  1.) 

(A  selection  of  this  information  is  provided  in  Appendix  A.) 

13.1  Enter  the  flag  number  =  "42"  upon  which  the  equilibrium  chemistry  read 
section  is  keyed  (if  you  use  the  kinetics,  you  must  include  all  species  from 
the  kinetics  section  in  the  equilibrium  section  in  order  to  load  the  thermo¬ 
chemical  arrays.) 

14.1  Enter  the  following  phrase  with  the  following  format  as  a  key  for 
chemistry  read  section: 

FORMAT ( A30,I3) :  "NUMBER  OF  REACTANT  ELEMENTS  =  " 

15.1  Enter  the  number  of  gas  phase  reactant  elements  considered  (-). 

16.1  Enter  the  following  phrase  with  the  following  format  as  a  key  for 
chemistry  read  section: 

FORMAT(A30):  "CONSIDERED  ELEMENTS  =  " 

17.1- 17.n  Enter  the  "n"  gas  phase  elements  considered  in  atomic  order 

FORMAT  (n(A4)) 

18.1  Enter  the  following  phrase  with  the  following  format  as  a  key  for 
chemistry  read  section: 

FORMAT(A30):  "REACTANTS" 

Enter  reactant  information  using  FORMAT  (5A4,4(A2,F3.0),A1) 

19.1- 19.5  Enter  the  reactant  species  name  (left  justified)  as  it  appears  in  the 

N-L  library. 

19.6-19.7  Enter  the  first  element  symbol  (left  justified)  and  the  number  of 
atoms  (right  justified). 

19.12-19.13  Enter  the  last  element  symbol  (left  justified)  and  the  number  of 
atoms  (right  justified).  Leave  any  extra  fields  blank. 

19.14  Enter  the  phase  as  "G" ,  "S",  or  "L"  for  gas,  solid,  or  liquid  state. 

19.15  Enter  Critical  Pressure  PCRIT  (Pa). 

19.16  Enter  Critical  Temperature  (K). 

19.17  Enter  SIGMA  for  Lennard  Jones  diffusion  A. 

19.18  Enter  Epsilon  over  K  EPSOK  for  Lennard  Jones  diffusion  (K). 

Continue  with  Card  19's  until  all  gaseous  reactants  have  been  entered. 
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20.1 


Enter  the  following  phrase  with  the  following  format  as  a  key  for 
chemistry  read  section: 

FORMAT (A30):  "PRODUCTS" 

21.1  Enter  the  considered  products  (gaseous,  liquid,  solidus)  in  the  same 
manner  as  the  reactants  with  the  same  format  until  all  products  are 
entered. 

Note:  To  match  internal  keyed  indices: 

The  9th  reactant  must  be  C(GR). 

The  10th  reactant  must  be  FE(A). 

The  11th  reactant  must  be  FE(C). 

The  15th  species  (#Reactants  +  #Products)  must  be  FE3C. 

Kinetics  data  input  follows 

C2  Enter  tire  following  phrase  with  the  following  format  as  a  key  for  kinetics 
read  section: 

FORMAT(A30):  "KINETICS  BEGIN  REACTION  CARDS" 

Contingent  upon  8.1  =  2. 

C.3.1  Enter  number  of  kinetics  reactant  species  (-). 

C.3.2  Enter  number  of  kinetic  reactions  (-). 

C.4.1  Enter  considered  reactant  species  FORMAT(4(5A4)). 

Enter  reactant  information  using  FORMAT(A10,7E10.4,I5). 

C.5.1  Enter  key  word  "REACTANTS". 

C.5.2  Enter  forward  rate  coefficient  CF  in  (cmA3/mole-s)  or  (cmA6/moleA2-s). 
C.5.3  Enter  forward  rate  temperature  exponent  ETAF  (-). 

C.5.4  Enter  forward  rate  activation  energy  EF/k  (K). 

C.5.2  Enter  backward  rate  coefficient  CB  in  (cmA3/mole-s)  or  (cmA6 
/moleA2-s). 

C.5.3  Enter  backward  rate  temperature  exponent  ETAB  (-). 

C.5.4  Enter  backward  rate  activation  energy  EB/k(K). 

If  backward  rates  are  entered  as  zero  then  they  are  not  considered,  if  they 
are  entered  as  "-1"  then  rates  are  determined  from  thermodynamic  data. 

C.5.5  Enter  third  body  flag  THIRDBODY  (-)  (1=T,  0=F). 
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C.5.6  Enter  surface  or  gas  phase  reaction  flag  ICOFU  (-) 

(0  NORMAL  1  =  SURFACE) 

Enter  reactant  and  product  data  using: 

FORMAT(5A4,F3.0,4(A2,F3.0),A1,4E10.4) 

C.6.1  Enter  reactant  of  first  reaction  using  NASA  Lewis  library  notation  (-). 
C.6.2  Enter  coefficient  of  first  atomic  in  reactant  (-). 

C.6.3  Enter  first  atomic  of  first  reactant  (-) . 

C.6.4  Enter  second  coefficient  of  second  atomic  in  reactant  (-)  blank  if  non¬ 
existent. 

C.6.5  Enter  second  atomic  of  first  reactant  (-)  blank  if  non-existent. 

C.6.6  Enter  third  coefficient  of  third  atomic  in  reactant  (-)  blank  if  non-existent. 

C.6.7  Enter  third  atomic  of  first  reactant  (-)  blank  if  non-existent. 

Enter  all  reactants  as  previous. 

C.7.1  Enter  key  word  "PRODUCTS". 

C.8.1-C.8.7  Enter  product  data  using  reactant  data  guide. 

CC.1.1  If  reaction  is  third  body,  enter  key  word  "EFFICIENCY". 

CC.1.2  Enter  efficiency  for  reaction  to  each  kinetic  reactant  species  in  the  order 
entered  in  C.4.1.  FORMAT(7(F10.0)). 

CC.2.1  If  reaction  is  not  a  third  body  reaction,  enter  key  word  "KEND". 

Continue  with  next  reaction  card  set  C.5-CC.2  until  reactions  are  completed. 


5.  ATEC  Carc^Variable  Definitions 


Card  1:  Initial  State  of  Bore  Gas 

TINTT  Initial  Gas  Temperature  (NOVA  run,  R) 

PINIT  Initial  Gas  Pressure  (NOVA  nm,  psi) 

UINIT  Initial  Gas  Velocity  (NOVA  run,  in/s) 

Card  2:  Flags,  Counters,  and  Grid  Layout 
CHROME  Thickness  of  chrome  layer  (in) 
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TMAXV 


Total  time  per  shot,  i.e.,  from  the  firing  of  one  shot  to  the  next 
shot(ms) 

XCON  Breech  velocity  augmentation  for  singularity  point  (-) 

XCHAMB  Chamber  length  (in) 

ZMAX  Total  barrel  length  (in)  =  (chamber  length  +  tube  length)  (in) 

NLOC  No.  of  axial  grid  points  (  <  100,  keep  larger  than  25) 

ZSTART  Start  of  distribution  of  100  axial  points  (in) 

ZEND  End  of  distribution  of  100  axial  points  (in) 

LP2D  Flag  for  2-D  LP  chamber  calculations  (-) 

This  function  is  NOT  IMPLEMENTED,  enter  0. 

LPOFFSET  Initial  piston  offset  for  2-D  LP  chamber  calculations  (in) 

This  function  is  NOT  IMPLEMENTED,  enter  0. 

IMC  Flag  to  cause  IMC  card  to  be  read  for  forced  convection 

boundary  condition  for  representing  integral  midwall  cooling 
(-) 

This  function  is  NOT  IMPLEMENTED,  enter  0. 

INSOBC  Flag  to  trip  insulated  outer  barrel  surface  (-) 

This  function  is  NOT  IMPLEMENTED,  enter  0. 

LBLOWBY  Flag  to  cause  blowby  multiplication  factors  card  to  be  read 
(-) 

This  function  is  NOT  IMPLEMENTED,  enter  0. 

AUGMENT  Convective  heat  transfer  multiplication  factor  1.0  or  user 
specified  value  (-) 

Conditional  Card  1:  Roughness  Data  (if  AUGMENT  is  negative) 

ESMOOTH  Smooth  deformation  limit  (Moody  Chart)  (m) 

EROUGH  Rough  deformation  limit  (Moody  Chart)  (m) 

Card  3:  Barrel  Geometry 

NDTZ  No.  of  axial  barrel  geometry  description  locations  (up  to  8) 

ZDB(l-NDTZ)  Axial  Location  (in) 

Card  4:  O.D. 

DBAR(l-NDTZ)  Barrel  O.D.  at  the  location  (in) 

Card  5:  Inner  Diameter 

DBOR(l-NDTZ)  Barrel  I.D.  at  the  location  (in) 
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Card  6:  Probe  Location  Information 

NHTZ  No.  of  probe  locations  (up  to  8) 

ZOUT(l-NHTZ)  Probe  AXIAL  LOCATION(s)  for  output  (in) 

Card  7:  Radial  Probe  Location 

ROUT (1-NHTZ)  Probe  RADIAL  LOCATION(s)  from  inside  surface  (in) 

Card  8:  Erosion  Material  (Base)  and  Numerical  Information 

ICHEM  1  equilibrium,  0  none,  2  finite  rate  kinetics  on  surface  and 

equilibrium  subsurface  (-) 

TFREEZ  Surface  routine  freeze  out  temperature  (K) 

RHOSTEEL  Density  of  the  wall  material  (kg/ m3) 

TMELT  Temperature  at  which  the  material  begins  to  melt  (K) 

CHTOUT  Outer  barrel  heat  transfer  coefficient  (W /  m2/  K) 

TOUT  Ambient  outer  air  temperature  (K) 

DETA  Computational  coordinate  spacing  (-)  usually  "1" 

HOF  Latent  heat  of  fusion  for  the  material  (J/kg) 

IAXI  "0"  if  the  system  is  planer  or  "1"  if  the  system  is 

axisymetric 

MTMIX  Surface  temperature 

(1)  Molecular  weighting  function 

(2)  Wall  temperature 

Conditional  Card  2:  Erosion  Material  (Coating)  and  Numerical  Information 
(if  chrome  >  0) 

NCOAT  Number  of  radial  points  in  coating  material  (-) 

NBASE  Number  of  radial  points  in  base  material  (-) 

TFREEZC  Coating  freeze-out  temperature  (K) 

RHOCOAT  Coating  density  (Kg/  m3) 

TMELTC  Melting  temperature  of  coating  (K) 

DETA  Computational  coordinate  spacing  (-)  usually  "1"  (-) 

HOFC  Latent  heat  of  fusion  for  coating  material  (J/kg) 

SCALC  Coating  Stress  Flag  (-) 
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Conditional  Card  3:  Stress  Calculation  Card  (if  SCALC  >=  1) 

YMC  Young's  Modulus  Coating  (Pa) 

CT EC  Coefficient  of  Thermal  Expansion  Coating  (1  /  K) 

PRC  Poisson's  Ratio  Coating  (-) 

USC  Ultimate  Strength  of  Coating  (Pa) 

YMB  Young's  Modulus  Base  (Pa) 

CTEB  Coefficient  of  Thermal  Expansion  Base  (1/K) 

PRB  Poisson's  Ratio  Base  (-) 

Conditional  Card  4:  Surface  Crack  Card  (if  SCALC  >  1) 

CW  Crack  Width  (m) 

CLIMIT  Crack  Distribution  Limit  for  Heat  Transfer  Augmentation  (m) 

Card  9:  Burst  Fire,  Numerical  Time  Step,  and  Print  Control 

NROUNDS  No.  of  rounds  to  be  simulated 

(NOT  IMPLEMENTED  IN  THIS  VERSION  =  1) 

NKOUNT1  Counter  for  output  print  intervals  (during  htemp.out  for  each 
round) 

NKOUNT2  Counter  for  output  print  intervals  (after  htemp.out  for  each 
round) 

DELTV  Numerical  time  step  (s) 

LASTSHOT  Last  round  of  the  immediately  previous  simulation 

(NOT  IMPLEMENTED  IN  THIS  VERSION  =  0) 

Card  10:  Output  File  Control  (0  =  false),  (1  =  true),  (2  =  >single  shot  =  burst) 


IUNIT 

Units  to  be  printed  0=K,  1=R 

ITEMP 

Temperature  0  =  as  computed,  1=  temperature  rise 

IOUT(l) 

Screen  output 

IOUT(2) 

Single  shot  probe  temperature 

(NHTZ  locations) 

IOUT(3) 

Single  shot  inside  surface  temperature 

(NHTZ  locations) 

IOUT(4) 

Single  shot  outside  surface  temperature 

(NHTZ  locations) 

IOUT(5) 

Burst  fire  probe  temperature 

(NHTZ  locations) 

IOUT(6) 

Burst  fire  inside  surface  temperature 

(NHTZ  locations) 

IOUT(7) 

Burst  fire  outside  surface  temperature 

(NHTZ  locations) 
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IOUT(8) 

Heat  transfer  coefficient 

(NHTZ  locations) 

IOUT(9) 

Effective  gas  temperature  history 

(NHTZ  locations) 

IOUT(IO) 

Single  shot  graphic  output  file 

(for  IBGRAF) 

IOUT(ll) 

Burst  fire  graphic  output  file 

(for  IBGRAF) 

IOUT(12) 

Boundary  condition  output  file  (P,T,U,  I.B.)  (for  IBGRAF) 

IOUT(13) 

Bore  surface  heat  flux  output  file 

(NHTZ  locations) 

IOUT(14) 

Bore  surface  regression  rate 

(NHTZ  locations) 

IOUT(15) 

Bore  surface  regression 

(NHTZ  locations) 

IOUT(16) 

Bore  surface  concentrations 

(NHTZ  locations) 

IOUT(17) 

Bore  core  flow  concentrations 

(NHTZ  locations) 

IOUT(18) 

Shear  stress  azmuthal  at  interface 

(NHTZ  locations) 

IOUT(19) 

Shear  stress  axial  at  interface 

(NHTZ  locations) 

(Note:  Screen  output  is  for  PC  applications  only  and  is  not  supported.) 
Card  11:  Barrel  Material  Properties 
IF4340  (1  =  true;  0  =  false) 

XKS  Conductivity  of  barrel  (W /m/K) 

(if  IF4340  =  1,  a  number  must  still  be  entered  here) 

ALFAS  Diffusivity  of  barrel  (m2/s) ) 

(if  IF4340  =  1,  a  number  must  still  be  entered  here) 

XKC  Conductivity  of  liner  (W/m/K) 

ALFAC  Diffusivity  of  liner  (m2/  s) 

Card  12:  Chemistry  Data  Flag  (-) 

AMAGIC  "42"  -  key  flag  number  for  the  chemistry  reading  section  (-) 

Card  13:  Chemistry  Data  Flag  (A30) 

TITLE:  Key  phrase  for  chemistry  read  section 

"NUMBER  OF  REACTANT  ELEMENTS  =  " 

Card  14:  Chemistry  Data  Card  (13) 

NEL:  Number  of  reactant  elements  (-) 
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Card  15:  Chemistry  Data  Flag  (A30) 

TITLE:  Key  phrase  for  chemistry  read  section 

"CONSIDERED  ELEMENTS  = " 


Card  16:  Chemistry  Data  Card  (10A4) 

CELEM2(NN) — Element  symbols 
Card  17:  Chemistry  Data  Flag  (A30) 

TITLE:  Key  phrase  for  chemistry  read  section 

"REACTANTS" 


Card  18  +  No.  of  Reactants:  Chemistry  Data  Cards  (5A4,4(A2,F3.0),A1) 


19.1-19.7 

CNAME(N,K): 

OSPECl(NJ) 

OSPEC2(NJ) 

ORFAZ(N) 

PCRIT(N) 

TCRIT(N) 

SIGMA(N) 

EPSOK(N) 


Enter  the  reactant  species  name  (left  justified)  as  it  appears 
in  the  N-L  library. 

Elemental  symbol  (left  justified) 

Number  of  elements  in  reactant 

Phase:  "G'V'L",  or  "S" 

Critical  Pressure  (Pa) 

Critical  Temperature  (K) 

Sigma  for  Lennard  Jones  diffusion  (A) 

Epsilon  over  K  EPSOK  for  Lennard  Jones  diffusion  (K) 


Continue  until  all  reactants  have  been  entered  #  =  MXG  AS 


Card  19:  Chemistry  Data  Flag  (A30) 
TITLE:  Key  for  chemistry  read  section 
"PRODUCTS" 


Card  20  +  No.  Products:  Chemistry  Data  Cards  (5A4,4(A2,F3.0),A1) 
20.1-20.7 


CNAME(N,K): 

OSPECl(NJ) 

OSPEC2(N,J) 

ORFAZ(N) 

PCRIT(N) 


Enter  the  reactant  species  name  (left  justified)  as  it  appears 
in  the  N-L  library. 

Elemental  symbol  (left  justified) 

Number  of  elements  in  reactant 

Phase:  "G",  "L",  or  "S" 

Critical  Pressure  (Pa) 
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TCRIT(N)  Critical  Temperature  (K) 

SIGMA(N)  Sigma  for  Lennard  Jones  diffusion  (A) 

EPSOK(N)  Epsilon  over  K  EPSOK  for  Lennard  Jones  diffusion  (K) 

Continue  until  all  products  have  been  entered  #  =  MXGAS 
Kinetics  data  input  follows: 

Conditional  Card  5:  Kinetics  Title  Card  (A40)  (if  ICHEM  =  2) 

TITLEK  Key  phrase  for  kinetics  read  section 


"KINETICS  BEGIN  REACTION  CARDS' 

Conditional  Card  6:  Kinetics  Data  Card  (315)  (if  ICHEM  =  2) 

KNSPEC  Number  of  species 

NRXN  Number  of  reactions 

Conditional  Card  7:  Kinetics  Data  Card  (4(5A4)  (if  ICHEM  =  2) 

CSPECK(NNJ)  Character  array  of  considered  species  names 

(Order  is  important) 

Conditional  Card  8:  Reactants  Title  &  Data  Card  (A40)  (if  ICHEM  =  2) 

TITLEK  Key  phrase  for  kinetics  read  section 

"CONSIDERED  SPECIES  FOR  KINETIC  REACTIONS" 
Conditional  Card  9:  Reactants  Name  Card  (4(5A4))  (if  ICHEM  =  2) 


CSPECK  Species  characters  using  NASA  Lewis  notation,  using  as 

many  lines  as  needed 

Conditional  Card  10:  Data  Card  (A10,6E10.4,F5.0,I5)  (if  ICHEM  =  2) 
TITLESMALL  Key  phrase  for  kinetics  read  section 
"REACTANTS" 


CF(N) 

ETAF(N) 

EF(N) 

CB(N) 

ETAB(N) 


Forward  rate  coefficient  (cmA3/mole-s)  or  (cmA6/ moleA2-s) 
Forward  rate  temperature  exponent  ETAF  (-) 

Forward  rate  activation  energy  EF/k  (K) 

Backward  rate  coefficient  (cmA3/mole-s)  or  (cmA6/ moleA2-s) 
"0"  -  not  considered,  "-1"  computed  from  thermodynamics 
Backward  rate  temperature  exponent  ETAB  (-) 

"0"  -  not  considered,  "-1"  computed  from  thermodynamics 
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EB(N), 


Backward  rate  activation  energy  EB/k(K) 

"0"  -  not  considered,  "-1"  computed  from  thermodynamics 
BODY3(N)  Third  body  flag  (1=T,  0=F) 

ICOFU(N)  Surface  or  gas  phase  reaction  flag  (0  NORMAL  1  =  SURFACE) 

Kinetic  reactant  and  product  data  cards  for  reaction  follow: 

Conditional  Card  11:  Reactants  Data  Card  (5A4,F3.0,4(A2,F3.0),A1,4E10.4) 

(if  ICHEM  =  2) 

CKNAME(N,M,K),K=1,5)  First  reactant  of  reaction  using  NASA  Lewis 

Notation 

CRCTT(N,M)  Reactant  atomic  coefficient 

OSPECKl(N,M,J)  Reactant  atomic  character 

OSPECK2(N/M,J)/J=l,4)  Corresponding  reactant  atomic  subscript 
ORFAZK(N,M)  Reactant  phase  "G","S","L" 

Continue  with  Card  11  until  all  reactants  in  reaction  are  entered. 

Conditional  Card  12:  Products  Flag  Card  (5A4,F3.0,4(A2,F3.0),A1) 

(if  ICHEM  =  2) 

CKNAME(N,M/K),K=1,5)  Key  phrase  for  kinetics  read  section 

"PRODUCTS" 

Conditional  Card  13:  Products  Data  Card  (5A4,F3.0,4(A2,F3.0),A1,4E10.4) 

(if  ICHEM  =  2) 

CKNAME(N,M/K)/K=1,5)  First  product  of  reaction  using  NASA  Lewis 

Notation 

CRCTT(N,M)  Product  Atomic  Coefficient 

OSPECKl(N,M,J)  Product  atomic  character 

OSPECK2(N,M,J),J=l,4)  Corresponding  reactant  atomic  subscript 

ORFAZK(N,M)  Product  phase  "G","S","L" 

Continue  with  Card  13  until  all  products  in  reaction  are  entered. 

Conditional  Card  14:  Efficiency  Flag  Card  (5A4,F3.0,4(A2,F3.0),A1) 

(if  BODY3  =  l) 

CKNAME(N/M,K),K=1,5)  Key  phrase  for  kinetics  read  section 

"EFFICIENCY" 
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Conditional  Card  15:  Efficiency  Data  Card  (7E10.4)  (if  BODY3=l) 

EFF(N,J)  Efficiency  of  each  species  CSPECK  (in  CSPECK  order!)  reacting 
with  each  other  species  (in  CSPECK  order!)  to  make  a  product. 

Conditional  Card  16:  End  of  Kinetic  Reaction  Card  (5A4,F3.0,4(A2,F3.0),A1) 

(if  BODY3  not  equal  to  1) 

KEND  Key  phrase  for  kinetics  read  section  after  each  reaction  if  it  is  not  a 

third  body 

"KEND" 

Continue  looping  on  conditional  Cards  10-16  until  all  reactions  are  entered. 
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Figure  1.  Description  of  post-mortem  erosion  pit. 
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Figure  2.  Analytical  surface  coating  description  including  stress. 
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Figure  3.  Code/file  connectivity. 
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Table  1.  ATEC  input/ output  file  table. 


"EK58.f"  file  numbering  table 

File 

No. 

File  name 

Purpose 

8 

htcoff.dat 

Heat  transfer  coefficient  data 

X 

9 

bprobe.dat 

Residual  probe  data  (burst) 

X 

10 

ssurfi.dat 

Single  (burst)  shot  bore 
temperature 

X 

11 

bsurfo.dat 

Residual  (burst)  outer  surface 
temperature 

X 

12 

LASTTEMP.DAT 

Restart  data  file 

X 

X 

15 

erod.in 

User  input 

X 

16 

erod.out 

Echo  of  user  input 

X 

17 

scratch 

— 

X  ! 

20 

sprobe.dat 

Single  (burst)  shot  probe 
temperature 

X 

21 

ssurfi.dat 

Single  (burst)  shot  bore 
temperature 

X 

22 

ssurfo.dat  i 

Single  (burst)  outer  surface 
temperature 

X 

23 

htemp.out 

IB  input 

X 

24 

LASTTEMP.DAT 

Restart  data  file 

X 

X 

25 

gastmp.dat 

Effective  gas  temperature  for 
exported  B.C. 

X 

27 

stemps.tab 

Single  (burst)  graphic  probe 
temperatures 

X 

28 

btemps.tab 

Burst,  graphic  probe  temperatures 

! 

X 

29 

sputib.tab 

Single  (burst)  graphic  P,T,u  output 
(tableaux) 

X 

30 

species.out 

Species  file  from  IBBLAKE 

X 

32 

library.in 

NASA  Lewis  Library 

X 

33 

— 

— 

— 

— 

34 

intp.out 

Interpolated  species  concentration 
output 

X 

35 


senergy.dat 


Bore  heat  flux  convection/ chemical 


Table  1.  ATEC  input/ output  file  table  (continued). 


"EK58.f"  file  numbering  table 

File 

No. 

File  name 

Purpose 

n 

36 

sdsdtdat 

Bore  surface  regression  rate 

X 

sregres.dat 

Bore  surface  regression 

X 

surfphi.dat 

Bore  surface  concentrations 

X 

39 

corephi.dat 

Bore  core  flow  concentrations 

X 

■ 

shears.dat 

Shear  stress  azmuthal  at  interface 

X 

!■ 

shearz.dat 

Shear  stress  axial  at  interface 

X 

fort.3 

Input  echo 

X 

testout 

Diagnostics  information 

X 
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Appendix  A.  Example  Input  Number  1 
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Example  Input  Number  1,  June  2001 
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....  Crack  Width  Cm) 

:  ....  Crack  Distribution  Limit  (m)  for  augmentation  of  HXF  in  crack 
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iunit  (0  for  K,  1  for  deg  R)  Units  flag 
_  itemp  (0  for  as  calculated,  1  for  temp  rise) 
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s:  The  variable  property  option  overrides  the  other  property  inputs  Alpa,Cp 
physical  property  inputs  are  in  units  of  W,m,s,K  (only  for  version  v32  and  above) 
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NUMBER  OF  REACTANT  ELEMENTS 
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Neither  the  critical  values  nor  the  Leonard  Jones  values  for  0  and  H  in 
substrate.  This  influences  neither  the  surface  state  nor  transport  calculation 
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EB/k (K) (0=>NOT, -1=>THERM0) 


; 1RST  REACTION  CO  +  0  <->  C02  3rd  order  rxn 
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;  FORMAT  (5A4 , F3 . 0 , 4 (A2 , F3 . 0 ) , A1 , 4E10 . 4 ) 
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EFFICIENCY 
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Appendix  B.  Example  Input  Number  2 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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ROUGHNESS  DATA  (only  if  heat  transfer  coef  fudge  factor  is  negative) 
esmooth  (m)  smooth  deformation  limit  (Moody  chart) 

:  erough  (m)  rough  deformation  limit  (Moody  chart) 
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OUTER  BARREL  HEAT  TRANSFER  COEF  (W/m2/K) 
....  AMBIENT  AIR  TEMP  (K)  (OUTSIDE) 
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186. D9  6.8D-6  0.35D0  900. D6  211. D9  12.6D-6  .296D0 
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4340  g.  362  curve  3  TPRC  Vol.10  0.057  cm2/s  @1203K  =  5.7e-6  m2/s 

CHROME  .  755  Inc.  and  Dewitt  0.207  cm2/s  @  600K  =  2 . 0795e-5m2/s 

CHROME  .  755  Inc.  and  Dewitt  0.148  cm2/s  @1000K  =  1 . 4828e-5m2/s 

Composite  Tzeng  Comp.  Manu.  Vol  4  no  3  0.0051  cm2/s  =  5.10e-7  m2/s 

Composite  Tzeng  Comp.  Manu.  Vol  4  no  3  0.1150  cm2/s  =  1.15e-5  m2/s 
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S02  s  10  2  G  7.88E6  430.8  4.112  335 

N2  N  2  G  3.39E6  126.0  3.798  71 

NO  N  10  1  G  6.48E6  180.0  3.492  116 

;  notes  the  1RST  solid  reactant  (9th  reactant)  must  be  C(GR)  solid  phase 
C (GR)  Cl  S  1.0  1.0  1.0  1 


the  10th  reactant  must  be  FE (A)  solid  phase 
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SIGMA,  EPSOK,  SOURCE  -  MARTIN  MILLER'S  BOOK  AND  HIRSHFIELDER 

TCRIT,  PCRIT,  SOURCE  -  VAN  WEHLAN  AND  SONNTAG 
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FE  (A)  1FE  1  S 
; note :  if  NOT  3rd  body  RXN,  KEND  is  required 
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EFFICIENCY 


;  7TH  REACTION  CO  +  02  <->  C02  +  O 
REACTANTS  5.06E13  0.0  31800. 
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Appendix  C.  Example  of  "species,  out"  File 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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Equilibrium  Compositions,  Rocket  Performance,  Incident  and  Reflected  Shocks,  and  Chapman-Jouget  Detonations."  NASA  SP- 
273,  NASA  Lewis,  Cleveland,  OH,  1971),  and  a  species  temporal  file,  "species. out",  from  a  coupled  interballistic  equilibrium 
thermodynamic  code.  The  ATEC  code  has  been  applied  to  the  25-mm  Bushmaster  cannon.  Navy  5-in  cannons,  the  120-mm 
M256  tank  cannon,  as  well  as  various  vented  fixtures,  and  has  performed  well  in  representing  the  physics  which  occurs  at  the 
surface. 
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